A Numerical Method for Parallel Particle Motions in Gyrokinetic Vlasov Simulations∗)
نویسندگان
چکیده
A semi-Lagrangian scheme is applied for the first time to computations of charged particle motions along magnetic field lines, to numerically solve the δf gyrokinetic equations in a flux tube geometry. This new solver adopted in the gyrokinetic Vlasov simulations has an advantage over the conventional Eulerian codes in calculating the parallel dynamics, because semi-Lagrangian schemes are free of the Courant-Friedrichs-Lewy (CFL) condition that restricts the time step size. A study of the accuracy of the parallel motion simulations reveals that numerical errors mainly stem from spatial (not temporal) discretization for realistic values of the grid spacing and time step, and it demonstrates the advantage of the semi-Lagrangian scheme. This novel numerical method is successfully applied to linear gyrokinetic simulations of the ion temperature gradient instability, where time steps larger than those restricted by the CFL condition can be employed.
منابع مشابه
An efficient numerical method for gyrokinetic δ f Vlasov simulations in helical plasmas
Three-dimensional configurations of helical plasmas bring numerical difficulties to gyrokinetic simulations of plasma turbulent transport in toroidal fusion devices. Watanabe et al. have investigated the ion-temperature-gradient (ITG)-driven turbulence in helical plasmas by means of the gyrokinetic Vlasov simulation code GVK [1]. Because of a large number of grid points along a magnetic field l...
متن کاملParticle-in-cell simulation with Vlasov ions and drift kinetic electrons
There are certain limitations in using gyrokinetic ions for simulations of turbulent transport in tokamak plasmas. Applications where Vlasov ions might be more appropriate include the electron temperature gradient driven turbulence, edge turbulence with steep density gradient, and magnetic reconnection in a weak guide field. In such cases the ion gyrokinetic model presently used in simulations ...
متن کاملTH/P8-30 Global Nonlinear Simulations of Ion and Electron Turbulence Usintg a Particle-In-Cell Approach
Ion-Temperature-Gradient and Trapped-Electron-Mode are two most likely candidates to explain the anomalous transport observed in Tokamaks. These turbulent phenomena are studied with the help of the gyrokinetic theory. These simulations require huge CPU power. The different numerical schemes all have some serious limitations. In the case of the Particle-in-Cell approach, the unavoidable statisti...
متن کاملComputation-Communication Overlap Techniques for Parallel Spectral Calculations in Gyrokinetic Vlasov Simulations
One of the important phenomena in magnetically-confined fusion plasma is plasma turbulence, which causes particle and heat transport and degrades plasma confinement. To address multi-scale turbulence including temporal and spatial scales of electrons and ions, we extend our gyrokinetic Vlasov simulation code GKV to run efficiently on peta-scale supercomputers. A key numerical technique is the p...
متن کاملGuiding-center simulations on curvilinear meshes
The purpose of this work is to design simulation tools for magnetised plasmas in the ITER project framework. The specific issue we consider is the simulation of turbulent transport in the core of a Tokamak plasma, for which a 5D gyrokinetic model is generally used, where the fast gyromotion of the particles in the strong magnetic field is averaged in order to remove the associated fast time-sca...
متن کامل